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摘 要 : 为 了 深入 发 掘 日 本 厚 朴 、 厚 朴 、 四 叶 厚 朴 叶 绿 体 基 因 组 差异 , 筛选 厚 朴 优良 性 状 候选 基因 ， 开 展 三 种 
厚 朴 的 分 子 遗 传 研究 ， 该 文 利用 Ilumina HiSeq 高 通 量 测序 平台 首次 对 日 本 厚 朴 叶绿体 进行 测序 、 组 装 ， 并 
与 已 有 的 厚 朴 、 四 叶 厚 朴 叶 绿 体 基因 组 共同 注释 ， 获 得 三 个 物种 叶绿体 基因 图 谱 ， 筛 选 出 三 个 基因 组 中 的 差 
异 基 因 ， 又 与 同 科 中 11 个 亲缘 物种 进行 叶绿体 基因 组 比 对 ， 构 建 NJ 遗传 树 。 结 果 表 明 : 日 本 厚 朴 叶绿体 基 
因 组 的 Clean Reads 为 19 791 019, Q30 73 91.3396, 组 装 后 基因 组 全 长 160 051 bp, GC 含量 为 39.2%, 含 tRNA 
37 个 ，rRNA 8 个 。 比 对 分 析 发 现 三 种 厚 村 具有 相似 的 IR LSC 和 SSC 结构 ， 以 及 GC 含量 和 tRNA 数量 
但 编码 基因 种 类 和 数量 、 内 含 子 和 外 显 子 的 数量 和 结构 等 存在 差异 。 日 本 厚 朴 的 功能 基因 数目 较 厚 朴 、 四 上 
厚 朴 分 别 多 6 个 和 4 个 ， 主 要 分 布 于 LSC 区 和 IR 区 ， 涉 及 核糖 体 大 亚 基 、 核 糖 体 小 亚 基 和 未 知 功能 基因 类 
群 。 系 统 发 育 分 析 结 果 进 一 步 显示 日 本 厚 朴 与 凹 叶 厚 朴 杂 缘 关系 较 近 ， 其 次 是 厚 朴 。 该 研究 表明 日 本 厚 朴 具 
有 更 丰富 的 叶绿体 基因 组 结构 、 组 成 和 变异 特征 ， 是 其 适应 高 纬度 地 区 弱 光 、 低 温 环 境 的 分 子 机 制 ， 这 为 厚 
朴 类 优良 品种 的 分 子 选 育 提供 有 力 的 指导 。 
关键 词 : 厚 朴 ， 凹 叶 厚 朴 ， 日 本 厚 朴 ， 叶 绿 体 基因 组 ， 系 统 发 育 树 
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Abstract: In order to investigate the good genes, cultivate the main superior cultivars and discover phylogenetic 
relationships of Magnolia officinalis, M. officinalis subsp. biloba and M. hypoleuca, we compared the differences 
among the cp (chloroplast) genomes of three Magnolia species and performed a phylogenetic tree of 14 species. 


Illumina HiSeq platform was used to sequence and assemble the cp genome of M. hypoleuca. Then the cp genomes 
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of three Magnolia species were annotated by online platform and performed with three Magnolia species cp gene 
cycles. Moreover, the cp genomes of other 11 Magnolia species were downloaded from the NCBI database and 
phylogenetic tree of 14 all species cp genomes was constructed based on NJ method. The results indicated that Clean 
Reads of M. hypoleuca were 19 791 019, and Q30 was 91.33%. The total length of cp genome of M. hypoleuca was 
160 051 bp, its GC content was 39.2%, including 37 tRNA and 8 rRNA. Besides, there were respectively six and 
four more functional gene numbers of M. hypoleuca than the other two Magnolia species, which indicated it had 
stronger viability. Compared with the cp genome structure of three Magnolia species, it showed that three Magnolia 
species had similar IR, LSC and SSC structures, GC content and tRNA number, but there were differences in the 
type and number of coding genes, the number and structure of introns and exons. Besides, the differential functional 
genes of three Magnolia species were mainly located in LSC region and IR region, involving large ribosomal 
subunits, small ribosomal subunits and unknown functional genes groups. Furthermore, based on NJ phylogenetic 
tree, M. hypoleuca was closely related to M. officinalis subsp. biloba, next to M. officinalis. In this study, M. 
hypoleuca had more abundant cp genome structure, composition and variation characteristics, which was the 
molecular mechanism of its adaptation to low light and low temperature environment in high latitude area. And it 
will also provide strong guidance for molecular breeding of excellent Magnolia varieties. 

Key words: Magnolia officinalis, Magnolia officinalis subsp. biloba, Magnolia hypoleuca, chloroplast genome, 
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厚 村 (Magnolia officinalis)、 四 叶 厚 机 CM. officinalis subsp. biloba) 及 日 本 厚 机 (M. hypoleuca) 均 为 木 
兰 科 药 、 材 两 用 经 济 林 阔 叶 树 种 ， 广 泛 分 布 于 恕 国 、 日 本 和 中 国 等 地 区 ， 有 具有 水 土 保 持 、 绿 化 、 美 化 环境 和 
药 用 等 多 种 应 用 价值 〈 彭 梦 婕 ，2020)。 目 前 三 种 厚 朴 的 相关 研究 主要 集中 在 化 学 成 分 、 临 床 应 用 及 资源 调 
E AZM, 2021) 等 方面 。 然 而 ,日 本 厚 朴 非 我 国 原 产 资源 ， 生 长 于 相对 低温 环境 中 ， 具 备 许多 优 于 另 
两 者 的 生物 学 性 状 ， 如 其 耐寒 能 力 相 比 于 其 他 木兰 属 植 物 更 强 CKwon & Oh, 2015)， 生 长 速度 和 成 熟 率 都 比 
同类 植物 更 快 ， 尤 其 是 在 生长 幼 期 ， 以 每 年 60~90 cm 的 速度 成 长 (Oguchi et al., 2017)， 是 解析 和 培育 厚 朴 
类 物种 优良 性 状 的 理想 材料 。 然 而 ， 日 本 厚 朴 的 分 子 遗传 学 研究 资料 荐 乏 ， 且 三 种 厚 朴 的 生长 速度 和 抗 寒 能 
力 差 异 的 形成 原因 缺乏 分 子 生 物 学 机 制 研究 ， 进 而 影响 了 厚 朴 类 优良 品种 的 选 育 。 

叶绿体 是 植物 细胞 中 必 不 可 少 的 细胞 器 ， 在 光合 作用 、 固 碳 以 及 氨基 酸 的 合成 中 起 到 重要 作用 ， 其 全 基 
因 组 包含 大 量 遗 传 信 息 ， 结 构 高 度 保守 ， 被 广泛 应 用 于 植物 分 子 进 化 及 系统 发 育 的 研究 ， 并 在 药 用 植物 的 遗 
传 转化 、 基 因 工 程 和 分 子 选 育 等 方面 发 挥 着 重要 作用 〈 赵 祺 等 ，2021)。 目 前 ， 厚 朴 和 四 叶 厚 朴 已 有 叶绿体 
基因 组 报道 , 而 具有 更 多 生物 特性 的 日 本 厚 朴 却 缺乏 叶绿体 基因 组 的 研究 。 鉴 于 此 , 本 研究 利用 Humina HiSeq 
高 通 量 测 序 技术 ， 对 日 本 厚 朴 进行 了 叶绿体 基因 组 测序 ， 与 既 有 的 厚 朴 、 止 叶 厚 朴 叶 绿 体 基 因 组 比 对 ， 算 选 
出 了 诸多 差异 的 核 苷 酸 序列 ， 破 译 三 种 厚 朴 之 间 的 亲缘 关系 ， 找 出 三 种 厚 朴 具有 更 强生 长 发 育 能 力 和 耐寒 性 
等 功能 的 相关 基因 ， 在 一 定 程 度 上 为 以 更 强 的 生长 速度 和 抗 寒 能 力 为 优良 性 状 的 目标 物种 提供 了 依据 。 
1 材料 与 方法 
1 材料 及 预 处 理 

日 本 厚 朴 新 鲜 、 幼 嫩 、 健 康 的 叶片 于 2019 年 10 月 10 日 采 自 四 川 省 成 都 市 温江 区 成 都 中 医药 大 学 药 用 


植物 园 (30°42'E，103°49'N)， 经 国家 中 药 种 质 资 源 库 专 家 高 继 海 副 教授 鉴定 所 有 样品 均 来 源 于 木兰 科 木 兰 
亚 属 日 本 厚 机 CM. hypoleuca)。 三 种 厚 机 的 命名 参考 于 中 国 植 物 物种 信息 数据 库 (http://db.kib.ac.cn/)。 叶 片 


二 


FEE 


7 


使 用 无 菌 水 擦拭 干将 
传统 文化 博物 馆 。 


叶绿体 提取 步 又 如 下 : orent Fr T RU 


1.25 mol.L NaCl, 0.25 mmol.L Vc, 1.5% PVP，pH3.6) (B, 4 


上 清 ; 加 入 常温 
DTT, 0.1% 牛 
存 备用 。 
1.2 叶绿体 基因 组 提取 和 测序 


清 


o ， 采 用 改良 的 CTAB (Matthes et al., 2020) 4 
EITA RAE TRÀ 
平台 Ilumina HiSeq PE150 进行 测序 


声波 机 械 打 断 ， 
后 使 有 


H 
市 


HH si ds 


， 迅 速冻 存 备 用 。 和 凭证 


8 接头 的 、 低 质量 的 Reads， 为 了 保证 
于 后 续 信息 分 析 。 数 据 过 滤 的 主 


标本 (馆藏 序列 号 为 ZYC190910) 保存 于 成 都 中 医药 大 学 中 医药 


研磨 后 悬 于 一 定量 A. 液 (50 mmol-L'! Tris, 25 mmol-L'! EDTA, 


出 修复 、 


ES 


置 ， 4 °C 下 2000 g 离心 10 min, SE Ej 


言 息 分 析 质 量 ， 对 Raw Reads 进行 过 滤 、 
要 步骤 如 下 : (1) 去 除 带 接头 的 Reads; 


层 纱 布 过 滤 ，4 °C 下 200 g 离心 20 min， 取 


缓冲 液 BC50 mmol-L'! Tris, 25 mmol-L'! EDTA, 1.25 mol-L'! NaCl, 0.25 mmol L'! Vc, 1 mmol L^! 
蛋白 BSA, pH 8.0)， 常 温 静 


， 叶 绿 体 沉 演 于 4 ?°C 保 


分 离 提取 DNA。DNA 经 检测 合格 后 ， 用 
溪 加 A、 连 接 测 序 接头 ， 进 行 PCR 扩 增 形成 测序 文库 ， 最 
字 。 测 序 得 到 的 原始 测序 序列 (Raw Reads)， 里 面 含有 
质 控 ， 得 到 Clean Reads， 用 
(2) 过 滤 N 含量 超过 10968] Reads; 


3'5 


, 


G) 去 除 质量 值 低 于 10 的 碱 基 超 过 50% 的 Reads。 对 过 滤 后 的 高 质量 数据 随机 抽取 2 000 条 Reads 数据 ， 
通过 Blast 软件 比 对 NT 库 检 测 样品 是 否 受到 污染 
1.3 叶绿体 基因 组 的 组 装 和 注释 
原始 序列 上 传 于 国家 生物 信息 中 心 数据 库 (序列 号 PRJCA004348 ) 。 利 用 Galaxy 在 线 平台 


(https://usegalaxy.org) 对 日 本 


时 村 叶绿体 基 


本 厚 朴 的 参考 基因 文件 。 再 
凹 叶 厚 朴 、 日 本 厚 朴 叶 绿 体 基 因 组 的 注释 。 
1.4 聚 类 分 析 
ERR. Pri ERR. XF NCBI 数 据 库 ， 


因 组 测序 结果 i 
NC_020316)〉 和 四 叶 厚 机 〈M. officinalis subsp. biloba, JN867581) 的 叶绿体 基因 
通过 CPGAVAS2 在 线 平台 


行 组 装 (Yan et al., 20150, 且 下 载 厚 机 CM. officinalis, 
组 FASTA 文件 ， 合 并 后 作为 


FS (http://www.herbalgenomics.org/cpgavas2) 完成 厚 朴 、 


下 载荷 花 玉兰 CM. grandiflora, JN867584)、 星 花 玉 兰 (Yulania 


stellata, NC_039941)、 天 女 木兰 (Oyama sieboldii, NC_041435)、 望 春玉 兰 (CY. biondii, KY085894)、 武 当 玉 


= (Y. sprengeri, JX280401)、 


yunnanensis, KF753638), 3E fÉ R75 CY. zenii, MH607378)、 紫 玉兰 CY liliiflora, X280397). $8 
chinense, NC. 0305040. JERE ZHK CL. tulipifera, NC. 008326) 
中 包括 2 REAREA 9 种 木兰 科 植 物 。 利 用 MEGAX 软件 ， 基 
内 的 14 种 木兰 科 植 物 的 系统 发 


2 结果 与 分 析 
2.1 叶绿体 基因 组 信息 


玉 灯 玉兰 CX 


通过 Ilumina HiSeq 高 通 量 测序 平台 


质量 的 Reads, 


兰 类 植物 的 最 大 3 
三 种 厚 朴 叶绿体 基 
厚 机) 和 88 146 bp < 


因 均 具 


测序 共 


获得 Clean Reads 19 791 019 条 ， 
bp， 符 合 目前 已 知 的 木兰 类 植物 叶绿体 基 
基因 组 (160 183 bp)， 而 


和 IRB) 分 别 为 26 566. 26 566 和 26 562 bp. 


E. E r 


中 GC 含量 最 高 
的 GC 值 最 小 ， 仅 为 34.296. 


的 区 域 都 为 下 区， 分别 为 43.2%、 
34.3% 和 34.3%。 
SSC WI IRA 区 ， 其 中 位 于 IRA 区 域 长 度 均 为 1 279 bp, FAM 


本 厚 朴 长 51 bp。 除 此 之 外 ， 


1279 bp. 
通过 基因 台 


结构 比较 ， 发 现 


育 树 ， 观 察 并 分 析 它 


因 组 大 小 范围 (159 429-160 183 bp)， 厚 朴 叶绿体 
[nm E 
有 典型 的 四 分 区 域 结构 ， 其 中 LSC 分 
本 厚 朴 )，SSC 分 别 为 18 843、18 832 和 18 771 bp，2 RE E48 
与 其 
本 厚 朴 特有 的 ycfl 


三 种 厚 村 叶绿体 的 rRNA 和 {tRNA 数量 相同 C 1)。 在 基因 数量 方面 ， 


云南 拟 单 性 木兰 CParakmeria 
3K Wk CLiriodendron 
木兰 科 植 物 的 叶绿体 基因 组 数据 ， 
邻接 法 (NJ 法 ) 构建 日 本 厚 朴 在 
门 之 间 的 亲缘 关系 (Yan et al., 2015). 


denudata ‘Lamp’, JN227740)、 


k 11 种 


N 


本 厚 朴 叶绿体 19 816 708. 条 原始 数据 ， 移 除 接头 和 低 
通过 在 线 组 装 发 现 日 本 厚 朴 全 长 160 051 
有 目前 已 知 木 


获得 
Q30 为 91.33%。 


的 叶绿体 基因 组 为 160 099 bp. 

) 别 为 88 210 bp《〈 厚 朴 )、88 145 bp (HN 
E 复 的 IR 区 (IRA 
—RUÉfIR. LSC 和 SSC 区 域 的 GC 值 存在 一 定 的 差异 ， 其 
43.1% 和 43.2%, LSC 区 次 之 ， 均 为 37.9%， 而 SSC KER 
由 物种 叶绿体 基因 组 相似 的 是 ， 三 种 厚 朴 的 ycf1 也 跨越 了 
[m EET SSC 区 域 长 度 为 4311 bp， 比 日 
x， 其 序列 长 度 分 别 为 29 bp. 


ZzET 


因 还 跨越 了 SSC 和 IRB 


pet 


个 单 拷贝 基 


厚 朴 相 比 厚 朴 多 2 


因 ， 而 两 者 又 较 


与 厚 朴 相 比 多 2 个 编码 区 ， 而 


本 厚 村 数量 


EE 最多， 比 前 两 者 分 别 多 8 个 和 10 个 。 这 揭示 了 厚 机 与 四 叶 


又 数量 , nr ER 


Ah 


本 厚 朴 分 别 少 了 4 个 和 6 个 。 在 编码 


量 方 


的 叶绿体 基因 组 可 能 存在 较 小 的 结构 和 功能 差异 ， 而 日 本 厚 朴 的 叶绿体 基因 组 结构 和 功能 可 能 较 其 他 两 者 更 
为 丰富 。 
de 1 三 种 厚 朴 叶绿体 基因 基本 信息 
Table 1 Basic chloroplast genes information of three Magnolia species 
单 拷贝 基因 总 数 
名 称 全 长 编码 区 
GC(%) Total number of single IRNA tRNA 
Name Length (bp) Coding region 
copy genes 
厚 朴 
160 183 39.2 126 82 8 37 
M. officinalis 
prim ER) 
M. officinalis subsp. 160 099 39.2 128 84 8 37 
biloba 
日 本 厚 朴 
160 051 39.2 132 92 8 37 


M. hypoleuca 


2.2 叶绿体 内 含 子 、 外 显 子 比较 


对 三 种 厚 朴 的 叶绿体 基因 


Z 


在 内 含 子 、 外 显 子 的 差异 而 3 个 
厚 朴 与 凹 叶 厚 朴 的 差异 主要 表现 在 11 个 ， 其 中 9 个 片段 的 差异 区 在 工 类 内 含 子 〈InI， 核 苷 酸 数量 差异 在 20 


以 内 )，2 个 在 工 类 内 含 子 〈InII[， 核 苷 酸 数量 差异 在 2 以 
EPII, 40 bp， 此 区 域 也 为 厚 朴 、 止 叶 厚 朴 叶 绿 体 基 
朴 存在 更 多 差异 ， 尤 其 是 trn1-GAU Æ 


能 基 


正常 表达 ， 可 


核糖 体 小 亚 基 基因 


[能 为 假 基因 


的 假 基因 


差异 可 为 厚 朴 类 杀 缘 关系 
K., IARTA IŽK, 


叶绿体 基因 组 条 


ZR 


clpP. ndhA. rpl2 等 基因 


注释 文件 进行 分 析 ， 发 现在 近 百 个 木兰 类 叶绿体 基 


基因 


因 的 1 类 内 含 子 比 前 两 者 多 920 ZER, HELE EF 
因 Crps22*, ycfI" fI ycf15^^) ( 表 2)， 这 些 结果 再 次 显示 厚 朴 与 止 叶 厚 朴 的 叶绿体 基 
小 ， 而 日 本 厚 朴 的 叶绿体 基因 


组 在 结构 和 功能 方面 较 其 他 两 者 更 为 丰富 。 


在 三 种 厚 朴 叶绿体 的 2 个 rps12 Æ 
。 植 物 叶 绿 体 基 因 
涉及 到 复杂 的 反 向 剪接 等 过 程 ， 本 研究 即 发 现 三 种 厚 朴 ; 
转化 在 物种 进化 过 程 中 起 到 重要 作用 
究 提 供 重 要 参考 。 三 种 厚 村 叶绿体 差异 最 大 的 ycf3 基 
II 类 外 显 子 ), 日 本 厚 村 的 I 类 外 显 子 增加 了 102 个 核 音 酸 。 此 外 ， 与 木兰 类 植物 的 
E43 trnaK-UUU、traL 等 相似 ,本 研究 针对 三 种 厚 朴 的 比较 分 析 , 发 现 trnI-GAU ,rps16、rpoC1、 
奶 位 点 率 较 高 ， 也 具备 开发 为 少数 木兰 植物 鉴定 条 形 码 的 潜力 。 
因 组 内 含 子 、 外 显 子 比较 


的 变异 位 点 率 、 信 


因 中 ， 


个 正常 含有 


表 2 三 种 厚 朴 叶绿体 基 


因 组 编码 基因 中 ，19 个 存 


的 内 含 子 、 外 显 子 数量 和 长 度 完全 一 致 ( 表 2)。 在 这 19 个 差异 片段 中 ， 


E 


内 )， 还 有 1 个 差异 片段 存在 开 类 外 显 子 Opl W 
办 组 差异 最 大 的 区 域 )。 相 较 于 前 两 个 厚 村 物种 ， 日 本 厚 


' 全 表达 的 功 
因 组 之 间 的 差异 较 


F 


73 个 外 显 子 ， 而 另 一 个 缺失 了 工 类 外 显 子 ， 无 法 
普遍 具有 保守 性 与 突变 性 〈 转 移 或 损失 等 ) 并 存 的 现象 ， 
rps12 的 几 个 外 显 子 为 反 向 排列 ， 且 为 非 顺 序 式 。 
(Liu et al., 2020)， 因 


其 转录 


此 rps12 核 苷 酸 序列 的 
K, 其 差异 区 域 数 达 4 个 (I 


Table 2 Comparison of introns and exons of three Magnolia species chloroplast genomes 


me 厚 朴 IDEE 本 厚 朴 
M. officinalis M. officinalis subsp. biloba M. hypoleuca 
dd EpI InI EpII Inl EpIl EPI InI EpII Ini  EpII Epl Inl Epi II EpIII 
trnK-UU 37 2498 35 37 2492 35 37 2493 35 
U 
rps16 42 824 246 42 824 246 40 823 221 
atpF 145 709 410 145 707 410 144 709 41 
432 740 1 614 432 734 1614 432 734 161 
rpoCl 4 
ycf3 124 733 232 729 151 124 734 232 727] 151 226 732 232 727 153 
trnL-UAA 35 491 50 35 491 50 35 491 50 
trnV-UAC 39 584 37 39 585 37 39 584 37 
rpl2^ 391 661 431 391 661 391 397 661 431 


rps12^ 114 536 232 26 114 526 232 26 114 536 232 26 
clpP 71 786 291 629 246 71 781 291 628 246 71 781 291 628 246 
petB 6 784 642 6 784 642 6 792 642 
petD’ 8 701 475 8 701 475 8 701 475 
rpll6 9 969 399 9 969 399 9 969 399 
ndhB^ 775 700 758 775 700 758 776 700 755 
trnl-GAU^ — 42 937 35 42 936 35 42 936 35 
trná-UGC 38 800 35 38 800 35 38 | 799 35 
b 
ndhA 553 1082 539 553 1102 539 552 1103 540 
trná-UGC 38 800 35 38 800 35 38 799 35 
trnl-GAU^ — 36 16 36 36 16 36 42 936 35 
ndhB^ 775 700 758 775 700 758 776 700 755 
rpl2* 391 661 431 391 661 431 397 661 431 
trnG-UC/ 
Mode 24 770 48 24 767 48 23 768 47 
WÈ: rps. 核糖 体 小 亚 基 ; rpo. RNA 聚合 酶 ; yef 开放 阅读 框 ; clp. 酷 蛋 白水 解 蛋白 酶 ;pet 多肽 ; ndh. NADH WAN: rpl. 核 
糖 体 大 亚 基 ; Ep. 外 显 子 ，In. AST; a. IRA X; b.IRB X; L.LSC 区 


Note: rps. Ribosomal protein subunit; rpo. RNA polymerase; ycf: Hypothetical chloroplast openreading frame; clp. Caseinolytic 


protease; pet. polypeptide; ndh. NADH dehydrogenase; rpl. Ribosomal protein large; Ep. Expressed region; In. Intron; a. IRA region; b. 
IRB region; L. LSC region. 


2.3 叶绿体 功能 基因 比较 


ME d 
CE 35, HEAT E 


基 


T 


在 叶绿体 四 大 类 基 
个 及 未 知 功能 基因 5A. PIER RS 
基因 


基因 


因 组 成 中 ， 


厚 村 含有 光合 作 有 


基因 表达 相关 基 


， 主 要 集中 于 


25k 
AEJ (Amborella), EXE 


wa (Nuphar), 


表达 相关 基 


昌 相 关 基 因 46 个 、 基 
因 少 2 个 (xps3 和 
ELA AR ARTE 
PV, CPszezólkowska et al. pen 功能 外 ， 


因 
rps15). ifj 


rE 


等 原始 被 子 植 


表达 相关 基因 6 
日 本 厚 朴 相 比 于 前 两 者 
EXED], BD rps3. rps15. ycfl. ycf15^^ 及 rpl22 
的 作用 还 存在 一 定 争 论 ， 如 ycf15 
物 中 是 无 功能 的 ， 甚 至 在 八 


、 其 他 基因 


69 个 


0 


角 属 CIllicium). ES (Acorus) fail (Ceratophyllum), p" EE (Ranunculus) 等 植物 的 进化 过 程 中 
已 经 完全 丢失 了 (Shi et al., 2013)， 而 木兰 属 (Magnolia) 和 胡椒 属 (Piper) 植物 基本 保留 了 这 类 基因 。 
综合 该 文 内 含 子 与 外 显 子 、 功 能 基因 的 结构 分 析 结果 , 三 种 厚 朴 叶绿体 基因 组 的 差异 主要 分 布 于 LSC 区 
和 IR 区 ， 涉 及 核糖 体 大 亚 基 、 核 糖 体 小 亚 基 和 未 知 功能 基因 类 和 群 ， 尤 其 是 诸多 未 知 功能 的 基因 ， 而 木兰 属 
中 这 类 疑似 非 功能 基因 的 结构 特征 、 变 异 原因 ， 以 及 与 适 生 环境 差异 的 关联 性 ， 有 待 进一步 深入 研究 。 
表 3 三 种 厚 朴 叶绿体 的 功能 基因 比较 
Table 3 Comparison of chloroplast functional genes of three Magnolia species 
基因 功能 基因 种 类 日 本 厚 朴 厚 朴 IDEE 
Gene function Gene type M. hypoleuca M. officinalis M. officinalis subsp. biloba 
ATP 合 酶 亚 基 
atpA, atpB, atpE, atpF, atpA, atpB, atpE, atpF, atpA, atpB, atpE, atpF, 
ATP synthase 
atpH, atpl atpH, atpl atpH, atpl 
subunit 
光合 系统 I 亚 基 
光合 psaA, psaB, psaC, psal, psaA, psaB, psaC, psal, psaA, psaB, psaC, psal, 
Photosynthetic 
作用 . psaJ psaJ psaJ 
system I subunit 
Photosynthesis 


光合 系统 II 亚 基 
Photosynthetic 


system II subunit 


psbA, psbB, psbC*, psbC*, 
psbD, psbE, psbF, psbl, 
psbJ, psbK, psbL, psbM, 
psbT, psbZ, psbH 


psbA, psbB, psbC*, psbD, 
psbE, psbF, psbl, psbJ, 
psbK, psbL, psbM, psbN, 
psbT, psbZ, psbH 


psbA, psbB, psbC*, psbD, 
psbE, psbF, psbl, psbJ, 
psbK, psbL, psbM, psbN, 
psbT, psbZ, psbH 


NADH 脱氧 酶 亚 基 


ndhA, ndhB^, ndhB^, 


ndhA, ndhB^, ndhB^, 


ndhA, ndhB^, ndhB^, ndhC, 


NADH ndhC, ndhD, ndhE, ndhF, | ndhC, ndhD, ndhE, ndhF, 
ndhD, ndhE, ndhF, ndhG, 
dehydrogenase ndhG, ndhH, ndhl, ndhJ, ndhG, ndhH, ndhl, ndhJ, 
ndhH, ndhl, ndhJ, ndhK 
subunit ndhK“, ndhK" ndhK 
细胞 色素 bf 复合 
体 petA, petB, petD*, petD', ~ petA, petB, petD", petG, petA, petB, petD', petG, 
Cytochrome b/f petG, petL, petN petL, petN petL, petN 
complex 
二 磷酸 核 酮 糖 羧 化 
KET. As 
基因 表达 
Ribulose rbcL rbcL rbcL 
Gene expression 
diphosphate 


carboxylase subunit 


核糖 体 大 亚 基 
Large ribosomal 


subunit 


rpl14, rpll6, rpl2", rpl2*, 
rpl20, rpl22, rpl23°, 
rpl23^, rpl32, rpl33, rpl36 


rpl14, rpl16, rpl2*, rpl2^, 
rpl20, rpi22, rpl23°, 
rpl23^, rpl32, rpl33, rpl36 


rpl14, rpl16, rpl2*, rpl2*, 
rpl20, rp122, rpl23^, rpI23^, 
rpl32, rpl33, rpl36 


核糖 体 小 亚 基 


Small ribosomal 


Subunit 


rps11, rps12°, rps12P, 
rps12*, rps14, rps15, 
rps16, rps18, rps19, rps2, 
rps3, rps4, rps7^, rps7", 
rpsá 


rps11, rps12^, rps12P, 
rps14, rps16, rps16, 
rps19, rps2, rps4, rps7", 
rps7,, rpsá 


rpsll, rps12*, rps12^, rps14, 
rps15, rpsl6, rps18, rps19, 
rps2, rps3, rps4, rps7', 
rps7,, rps& 


依赖 于 DNA 的 
RNA 聚合 酶 亚 基 
DNA-dependent 


rpoA, rpoB, rpoCl, 


rpoA, rpoB, rpoC1, 


rpoA, rpoB, rpoC1, rpoC2 


rpoC2*, rpoC2* rpoC2 
RNA polymerase 
subunit 
rrn16^, rrn23", rrn23°, rrn5?, rrn23^, rrn23^, 
" rrn16*, rrn23", rrn4.5", ; 5 ; 
核糖 体 RNA $ rrn4. 5", rrn5*, rrn5", rrnl16^, rrn4.5^, rrn16^, 
rrn5*, rrn5", rrn23", ' , i 
Ribosomal RNAs " à rrn23^, rrn23^, rrn16^, rrn23", rrn23", rrn4.5*, 
rrn16^, rrn4.5 A 
rrn4.5 rrn5*, 
乙酰 辅酶 A ERN 
亚 基 
accD accD accD 
Acetyl-CoA 
carboxylase subunit 
C 型 细胞 色素 合成 
ccsA ccsA ccsA 
C-type cytochrome 
synthase 
iata A A A 
cem cem cem 
其 他 基因 Membrane protein 
Other genes 和 蛋白酶 
clpP clpP clpP 
Protease 
翻译 起 始 因子 
Translation infA infA infA 
initiation factor 
成 熟 酶 matK matK matK 


202108.00056v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


Mature enzyme 


未 知 功能 基因 


Unknown 


functional gene 


保守 的 开放 阅读 杠 


Conservative open 


reading frame 


ycf1*^, ycfl A ycf2*, 
ycf2^, ycf3, ycf4, ycf15*, 
ycf15" 


ycf4 


yef1 **, ycf2", ycf2", ycf3 


yfl, ycf2", yef2", ycf3 
ycf4 


iE: psa. 光合 器 ; acce. 乙酰 
mat. 成 熟 酶 基因 |; 


甫 酶 A 羧 化 酶 ; 


rbe. 二 磷酸 核 酮 糖 羧 化 酶 ; 


ccs. C 型 细胞 色素 合成 卉 


Xl; cem. HRAJE 


a. IRA X; b.IRB X; s-a. 横 跨 SSR XAI IRA X; s-b. TEES SSR 


X fll IRB 


义 。 


Note: psa. Photosynthetic apparatus; acc. Acetyl-coa carboxylase; rbe. Ribulose bisphosphate carboxylase; ccs. TypeC cytochrome 


synthesis gene; cem. Chloroplast envelope membrane protein; mat. Mature enzyme gene; a. 1 IRA region; b. IRB region; s-a. Across the 


SSR and IRA regions; s-b. Across the SSR and IRB regions. 


2. 4 部 分 木兰 科 植物 的 亲缘 关系 


类 群 中 ， 


厚 朴 。 木 兰 属 中 另外 一 个 次 级 类 条 


由 图 1 可知， 在 14 种 近 缘 物种 叶绿体 基 


和 包括 云南 拟 贞 


性 木兰 、 玉 灯 玉 兰 、 紫 玉兰 、 


因 组 构建 的 NJ 发 育 树 〈 节 点 支持 率 均 在 75% 以 上 ) P, ER 
属 和 其 他 属 是 木兰 科 中 的 两 个 独立 的 谱系 ， 这 与 之 前 报道 的 结果 一 致 〈 陈 凯 ，2019)。 在 木兰 属 的 2 个 次 级 
三 种 厚 朴 、 天 女 木 兰 、 和 荷花 玉兰 聚 类 在 一 起 ， 其 中 日 本 厚 朴 与 凹 叶 厚 朴 植物 关系 最 为 密切 ， 


其 次 是 


武当 玉兰 、 望 春玉 兰 、 星 华 玉 


兰 、 宝 华 玉 兰 ， 它 们 与 三 种 厚 朴 亲缘 关系 由 近 及 远 。 


3 讨论 与 结论 


一 ) 


宝 华 玉兰 Wuüania zenii MH607378 


E 1G X Y stellataNC 039941 
BEEŽ Y biondiiKY085894 

武当 玉兰 sprengeri JX280401 
AE Y liliiflora X280397 

玉 灯 玉兰 了 denudata ‘Lamp’ JN227740 


云南 拟 单 性 木兰 Parakmeria yunnanensis KF753638 


荷花 玉兰 Magnolia grandiflora JN867584 
天 女 木兰 Oyama sieboldii NC 041435 

" Magnolia officinalis NC 020316 
凹 叶 厚 朴 M. officinalis subsp. biloba JN867581 


日 本 厚 朴 M. hypoleuca 


13 x tM Liriodendron chinense NC 030504 
JIŽ MEL. tulipifera NC 008326 


图 1 部 分 木兰 科 物 种 系统 进化 树 


Fig.1 Phylogenetic tree of partial Magnoliaceae plants 


EN 


究 使 用 全 基因 组 测序 来 组 装 和 比较 三 种 厚 朴 的 叶绿体 基因 
性 状 差 异 的 形成 原因 。 低 温 环 ] 


R, CO; 同化 作用 和 光 系 统 II 活 性 ， 从 而 影 


psbC 和 psbD 是 光合 系统 [I 的 关键 基 


响 植物 的 了 


组 ， 以 此 来 


测 三 种 厚 朴 之 间 多 种 生物 学 
啼 下 ， 植 物 最 先 受到 抑制 的 生理 代谢 过 程 是 光合 作用 ， 低 温 能 显著 降低 光合 效 
E 常 生长 发 育 ( 王 璐 等 ，2020)。 通 过 三 种 厚 村 叶绿体 的 
功能 基因 比较 ， 发 现 日 本 厚 机 psbC 基因 在 数量 上 较 另 两 者 增加 了 一 个 ， 这 可 能 缘 于 psbC 基因 的 自我 复制 。 
因 ， 例 如 小 麦 成 熟 叶 绿 体 psbC 和 蛋白 的 合成 发 生 于 含 psbD 序列 的 转录 本 
E, ÆR psbD-psbC 复合 物 ， 它 们 的 转录 水 平 受到 光 诱 导 ， 通 过 增强 光合 系统 I 避 
高 叶绿体 的 合成 速度 (Gamble et al., 1988; Gamble & Mullet, 1989)， 从 而 提高 叶绿体 的 光合 作 /) 


d= 


5 基 的 合成 和 维持 能 力 ， 提 


强度 ， 进 而 


加 快 植物 的 生长 发 育 。 此 外 ，psbD-psbC 基因 协同 转录 的 mRNAs 可 以 翻译 产生 D2 和 CP43 蛋白 ， 与 反应 中 


心 D1、CP47、 放 氧 复合 体 蛋 白 及 捕 光 复合 体 I SERE 


， 共 同 参与 叶绿体 光合 系统 I 的 光合 电子 传递 ， 在 强 


光 条 件 下 维持 光合 系统 工 的 功能 ( 庄 烛 扬 ，2020)。psbD-psbC 复合 体 直 接 参与 光合 系统 I[[ 亚 基 的 形成 ， 在 生 


长 旺盛 和 存在 颜色 差异 的 叶片 中 发 挥 着 保护 光合 系统 、 减 少 强 光 伤害 的 作用 
矢 竹 叶 等 光合 系统 I 中 psbD 基因 的 研究 中 也 得 到 证 实 〈 许 冰 清 等 , 2015 )。 因 
基因 表明 其 具备 合成 更 多 psbD-psbC 复合 体 与 D2 蛋白 的 潜力 ， 这 可 能 是 日 本 厚 朴 生长 速度 更 快 、 适 应 较 北 


纬度 地 区 低温 和 强 光 胁迫 环境 的 原因 之 一 


>H 


在 高 变异 位 点 率 的 功能 基因 中 , rpl22 为 日 本 厚 朴 中 特有 , 其 位 于 LSC 区 , 属于 基因 表达 相关 功能 基因 ， 


(Adachi et al., 2012)， 这 在 花 叶 
此 ,日 本 厚 朴 叶绿体 增加 的 psbC 


可 用 于 物种 鉴定 〈Feng et al., 2019)。 除 此 外 ， 本 研究 发 现 三 种 厚 朴 的 差异 基因 主要 分 布 于 rp122 所 属 的 基因 
IRR, CMER rp123 操纵 子 的 大 转录 单元 ， 翻 译 为 核糖 体 大 蛋白 亚 基 。 该 研究 发 现 ， 此 操纵 子 在 三 种 厚 朴 
中 存在 显著 的 结构 差异 : rpl23-rp12-rps19-rp122. C HA EE RO -rps3 CIL PE EE $h, El ZR JE kh) 
-rp116-rp114-rps8-infA-rpl36-rps11-rpoA , 3X —^ VPRBURZROS AR upAd, MHE HRERS EEk 
过 程 中 ， 该 操纵 子 rps19-rp116 Z WRA EZ TERRAS EOERUCC rp13. rp122 基因 填充 ， 同 时 造成 操纵 子 中 多 个 
功能 基因 的 内 含 子 、 外 显 子 核 苷 酸 发 生 了 增 减 《减少 为 主 )， 功 能 蛋白 的 数量 和 活性 也 发 生 了 改变 。 植 物 叶 
绿 体 除 了 光合 作用 外 ， 还 部 分 参与 氮 基 酸 、 核 音 酸 、 脂 类 和 淀粉 等 各 类 成 分 的 生物 合成 ， 支 撑 起 该 植物 的 生 
物 学 性 状 CNamgung et al., 2021), 日 本 厚 村 具备 更 丰富 的 核糖 体 基 因 徐 ,这 无 疑 形成 了 其 更 丰富 的 生物 学 特 
性 。 


tRNA 作为 核酸 信息 水 平和 蛋白质 功能 水 平 的 适配器 ， 在 蛋白 质 翻译 中 起 着 核心 作用 ， 其 结构 修饰 影响 
植物 体 的 温度 适应 能 力 (Lorenz et al., 2017)。 本 研究 发 现 日 本 厚 村 tRNA-ALA 比 其 他 两 种 厚 朴 多 3 个 ， 具 有 
快速 转运 和 积累 丙 氨 酸 的 潜力 ， 而 游离 态 的 丙 氨 酸 能 抵抗 多 种 外 界 环境 刺激 CMustroph et al., 2014)， 如 寒冷 
因子 。 除 了 数量 因素 外 ， 三 种 厚 朴 叶绿体 中 tRNA 对 应 的 氨基 酸 种 类 也 不 同 。 日 本 厚 朴 为 本 属 植 物 中 分 布 较 
为 北端 〈 如 千岛 群岛 ) 的 物种 ， 其 生境 中 弱 光 和 低温 信号 交叉 影响 着 植物 的 生长 发 育 ， 在 这 种 光 温 胁迫 条 件 
F, tRNA 可 能 一 方面 积累 更 多 的 游离 丙 氮 酸 ， 另 一 方面 降低 正常 蛋白 质 翻 译 速 率 ， 将 正常 蛋白 质 合成 转向 
应 激 反应 蛋白 质 的 合成 〈 汉 德江 等 ，2002)， 最 终 使 得 植株 具备 更 强 的 耐寒 性 ， 这 又 为 厚 朴 优良 品种 的 分 子 
筛选 和 培育 提出 了 参考 。 
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